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p • Abstract 

We present the complete bosonic contributions to the effective weak mixing angle, 

sin 2 at the two-loop level in the electroweak interactions. We find their size to 

be about three times smaller than inferred from simple estimates from lower orders. 
In particular, for a Higgs boson mass, Mh, of 100 GeV they amount to 4 x 10 -6 , 
and drop down by about an order of magnitude for Mh = 200 GeV. We estimate 



the intrinsic error of the theory prediction of sin^ 0^q^ to be 4.7 x 10 



1 Introduction 



While no clear experimental evidence for the Higgs boson has been found so 
far, even today the Standard Model is accurately tested and the Higgs boson 
mass Mh strongly constrained through precision measurements of Z- and W- 
boson properties. One of the most important observables in this context is the 

effective leptonic weak mixing angle sin 2 . It can be defined through the 
vertex form factors for the vector and axial-vector interactions between the Z 
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boson and fermions /: 

eS 4 ^ \9a(M}))) 



with the g Vj A the effective couplings in the vertex i l^{gv — 9Alb)l Z^. 

The experimental value for sin 2 9^^ is derived from various asymmetries mea- 
sured at the Z resonance pole. The precision of the current experimental value 

sin 2 9 g?* = 0.23153 ± 0.00016 [1] could be improved by about one order of 
magnitude by a future linear collider experiment [2] . Since the radiative correc- 
tions to sin 2 depend sensitively on the value of Mh, the high experimental 
precision allows to put strong constraints on the Higgs boson mass when the 
Standard Model is assumed to be valid. Thus a lot of effort has been put 

into accurate theoretical calculations for sin 2 9^^. While one-loop corrections 
and two- and three-loop QCD corrections have been known for several years 
[3,4,5], only recently the fermionic two-loop corrections, i.e. the two-loop con- 
tributions with at least one closed fermion loop, were computed in [6,7,8] and 
confirmed in [9]. In addition, leading three-loop effects of order 0(a 3 ) and 
0(a 2 a s ) for large values of the top quark mass m t have been calculated [10], 
as well as the behavior of the full 0(a 3 ) corrections for large M H [11]. Finally, 
the precision of the QCD corrections to the universal part, the Ap parameter, 
has been pushed to the four-loop level [12,13]. 

However, for the remaining bosonic two-loop corrections, only a partial result 
for the M#-dependent diagrams is available so far [14]. The goal of this work 
is to finalize the calculation of the 0(a 2 ) two-loop corrections by giving a 
complete result for the bosonic two-loop contributions. 

At tree-level, the effective weak mixing angle is identical to the on-shell weak 
mixing angle sin 2 6^ = 1 — M w jM\. The effect of higher-order corrections to 
the Zll vertex can be summarized in the quantity An, 

=(l-g)(l + A«), (2) 



where for the purpose of this work, it is understood that M w and M z are 

defined in the on-shell scheme. The most precise result for sin 2 9^^ is ob- 
tained when using the Fermi constant G M instead of Mw as input. Then the 

calculation of sin 2 9^^ as a function of G M involves also the computation of 
the radiative corrections to the relation between and Mw- This has been 
carried out with complete electroweak two-loop corrections in [15,16,17]. In 
this letter, the remaining bosonic two-loop corrections to the form factor An 



2 



are presented. 



2 Outline of the calculation 



Any higher order calculation consists of two parts: the computation of the bare 
diagrams and the determination of the renormalization constants. The latter 
has been discussed at length in connection to two-loop electroweak precision 
observables in [15,16] (see also [18]). We are, therefore, left with the calculation 
of the bare diagrams, which in our case, are massive two-loop three-point 
functions with two massless and one massive external leg. 

Just as in the case of the fermionic corrections [6] , there are three mass scales 
in the problem, with the difference that there is no dependence on the top 
quark, but on the Higgs boson mass. This is, however, an important difference, 
because contrary to m t , M H is not a fixed parameter and can assume a broad 
range of values. From the many possible strategies that one might apply, we 
chose to expand in the various parameters in order to obtain a result expressed 
through single scale integrals, which are in fact just numbers to be determined 
in a final step. 

In a first step, we apply an expansion in the difference of the masses of the 
W and Z bosons, where the expansion parameter is just s^. Since there are 
diagrams where there is a threshold when Mw = Mz, the appearance of 
divergences at higher orders in the expansion is inevitable. In this case, we 
apply the method of expansions by regions, see [19]. The two regions that 
contribute to the result come from the ultras oft momenta, k lj2 ~ s^Mz, and 
hard momenta, ~ Mz- The new integrals that appear from this procedure 
are presented in Ref . [20] , whereas the reduction to the set of master integrals 
proceeds with Integration-By-Parts identities [21] solved with the Laporta 
algorithm [22] as implemented in the IdSolver library [23] . 

The Higgs boson is treated in two regimes. For low masses we expand in the 
mass difference between Mh and Mz, with the expansion parameter defined 
to be 



9 Ml 

4 = 1-4, (3) 



where this time, no thresholds are encountered. To guarantee a reasonable 
precision, we compute six terms in the combined expansion in s\^ and s 2 H . For 
the second regime, which is the region where Mh 3> Mz, we apply the large 
mass expansion, Ref. [19]. 
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The resulting single scale master integrals are treated with various methods, 
usually with two or three different ones for test purposes. Most can be obtained 
with numerical integration, using dispersion relations (see Method II in [8]). 
For diagrams of simpler topologies we use differential equations [24,25] and 
large mass expansions, whereas for more complicated ones we used the MB pack- 
age [26] implementing Mellin-Barnes methods [27,28] (see also [29]). Whenever 
possible we performed cross checks with sector decomposition, Ref. [30]. 

A final, algebraic check of all the procedures is the cancellation of the depen- 
dence on the gauge parameter, which we verified for the first orders of the 
expansion. 



3 Results 

The different electroweak contributions to Ak are shown in Tab. 1, under the 
assumption that the W boson mass is fixed at its experimental value given in 
Tab. 2 together with the remaining input parameters. Note that A«( a2 ^ erm 
denotes the two-loop contribution of the fermionic diagrams known from [6,9], 
whereas Ak^^^ 008 is our new result, namely the two-loop contribution of 
the bosonic diagrams. The complete prediction, Ak, contains additionally 
other known corrections. In particular, the following have been taken into 
account (see also [6]): one-loop electroweak corrections, QCD corrections to 
the one-loop prediction at the two- [4] and three-loop level [5], 0(a 2 a s mj) 
and <D(a 3 mt) corrections to Ap [10], as well as leading reducible effects at 
0(a 2 a s ) and 0(a 3 ). The exact M H dependence of the two-loop contributions 

to sin 2 ^^, obtained from Ak by rescaling with 1 — M^/Mf, can also be 
read off from Fig. 1, which makes it even more apparent that both fermionic 
and bosonic corrections are of the same order for low to moderate Higgs boson 
masses. 
Table 1 

Higgs boson mass dependence of Ak evaluated with a fixed W boson mass as in 
Tab. 2. The normalization factor is 10~ 4 . 



M H [GeV] 




AK (a 2 ),ferm 


AK (a 2 ),b0S 


Ak 


100 


413.325 


1.07 


-0.74 


372.93 


200 


394.023 


-0.32 


-0.47 


353.20 


600 


354.060 


-2.89 


0.17 


313.13 


1000 


333.159 


-2.61 


1.11 


295.11 



In order to partially cancel large perturbative effects and lower the sensitivity 
to input parameters, the W boson mass is customarily replaced by its value 
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Fig. 1. Two-loop electroweak contributions to the effective weak mixing angle, with 
a fixed W boson mass as in Tab. 2. 

Table 2 

Input parameters, taken from [1,31]. 



input parameter 


value 


M w 


80.404 ± 0.0030 GeV 


M z 


91.1876 ± 0.0021 GeV 


Tz 


2.4952 GeV 


m t 


172.5 ± 2.3 GeV 


m b 


4.85 GeV 


Aa(Mf) 


0.05907 ± 0.00036 


a s (M z ) 


0.119 ±0.002 




1.16637 x 10~ 5 GeV" 2 



determined from [i decay, or equivalently from the Fermi constant, G^. The 
transition is made possible by the knowledge of all relevant corrections to 
muon decay as discussed in [17]. 

In view of the reparametrization, let us define the size of the complete bosonic 

corrections as the difference between the complete prediction of sin 2 6^ * and 
the same prediction, where pure two-loop electroweak bosonic diagrams have 
been omitted both in Mw and in Ak. A rough estimate of the effect can 
be read off Tab. 1. For example, for Mu = 100 GeV and the current input 

parameters, where sin 2 9w = 0.2225, the contribution to sin 2 g? * from Ak is 
—0.16 x 10 -4 , whereas using our previous results, [15,16,17], we find that the 
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contribution from Mw amounts to 0.2 x 10~ 4 . The large cancellation in the 
sum gives just 0.04 x 10~ 4 . A similar cancellation occurs for other values of 
M H over a wide range, as illustrated in Tab. 3, which also gives the complete 
prediction. It is important to note that, contrary to the fermionic corrections, 
which strongly depend on the value of the top quark mass, our result for the 
bosonic corrections is stable within the input parameter uncertainties and can 
be simply added to our fitting formula from [6], although the effect is clearly 
negligible. 

Table 3 

Higgs boson mass dependence of Mw and sin 2 6^ with as input parame- 
ter. Quantities with the superscript [ferm] do not contain the two-loop electroweak 
bosonic corrections, whereas those with the superscript [ferm+bos] do contain them. 
A sin 6> e u is the shift induced by the bosonic corrections as described in the text. 

The normalization factor for sin 2 0^3?* is 10~ 4 , and Mh and Mw are given in GeV. 



M H 


[fermj 

1V1 W 


[ferm+bos] 
M W 


sm e eS 


• 2 Jferm+bosj 
sm tf eff 


Asin 2 eff 


100 


80.3694 


80.3684 


0.231434 


0.231438 


0.04 


200 


80.3276 


80.3270 


0.231769 


0.231769 


0.00 


600 


80.2491 


80.2490 


0.232322 


0.232327 


0.05 


1000 


80.2134 


80.2141 


0.232563 


0.232574 


0.12 



As a partial check of our calculation we have also compared the Higgs boson 
mass dependence of our result with the one published in [14]. Upon using the 
same input parameters, we have found excellent agreement as shown in Tab. 4. 

Table 4 

Comparison of the Higgs boson mass dependence of An with [14] (with Mw = 
80.4260 GeV and nit = 178 GeV taken from that work). The subscript "sub" denotes 
subtraction at Mh = 100 GeV. The normalization factor is 10~ 4 . 



M H [GeV] 


A/ > 2 ),bos 
sub 




100 








200 


0.266 


0.265 


600 


0.914 


0.914 


1000 


1.849 


1.849 



4 Discussion 

Our calculation shows that the last piece of the two-loop electroweak correc- 
tions to the effective weak mixing angle, the one coming from purely bosonic 
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diagrams, gives a very small contribution. In fact, being of the order of a 
few times 10 -6 , it is below the anticipated precision of the linear collider, not 
even to mention the current experimental accuracy. This strengthens the va- 
lidity of our fitting formula [6]. Furthermore, since the recent calculation of 
the (9(G M m 2 a|) corrections to the rho parameter [12,13], has also given a very 
small contribution, the results of [6] implemented in ZFITTER [32] still provide 

a reliable prediction for sin 2 O^g*" under consideration of all new results. 

We estimate the error from unkown higher order corrections on sin 2 O^q^ as 
described in [6], including contributions for the next missing loop orders, i.e. 
C(a 3 ), 0(a 2 a s ), 0{aaf). 1 We find a total theoretical error of 4.7 x 10 , 
which should be taken as a very conservative error estimate, which makes 
it however necessary to determine further missing corrections if we want to 
obtain a result at the level needed by a future linear collider. 
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